Direct Mott Insulator-to-Superfluid Transition in the Presence of Disorder 
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We introduce a new renormalization group theory to examine the quantum phase transitions upon 
exiting the insulating phase of a disordered, strongly interacting boson system. For weak disorder 
we find a direct transition from this Mott insulator to the Superfluid phase. In d > 4 a finite region 
around the particle-hole symmetric point supports this direct transition, whereas for 2 < d < 4 
perturbative arguments suggest that the direct transition survives only precisely at commensurate 
filling. For strong disorder the renormalization trajectories pass next to two fixed points, describing 
a pair of distinct transitions; first from the Mott insulator to the Bose glass, and then from the 
Bose glass to the Superfluid. The latter fixed point possesses statistical particle-hole symmetry and 
a dynamical exponent z, equal to the dimension d. 
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The field of quantum phase transitions in disordered 
systems experienced remarkable growth recently. Its 
experimental relevance comprises the superconductor- 
insulator transition in thin films [0, "^He in disordered 
media 0, as well as Quantum Hall systems, and many 
types of disordered magnets Q . Quantum criticality may 
lie at the heart of the anomalous normal state properties 
of high Tc superconductors, and it can describe the char- 
acteristics of vortex systems with correlated disorder M . 
The intensely studied localization transition of interact- 
ing, disordered bosons has served as a useful paradigm for 
quantum phase transitions. Most of the originally pro- 
posed theoretical picture [p| was subsequently confirmed 
numerically H. However a key issue remains controver- 
sial. Analytic theories suggest, that the localized, "Bose 
Glass" phase intervenes everywhere between the interac- 
tion driven (Mott) insulator and the superfluid [|| . How- 
ever, none of the numerical studies of the problem [0j|] 
found the Bose glass in 2d at commensurate densities. 
Experimental tests in Josephson junction arrays also an- 
alyzed the critical behaviour [|lO[ . In these systems charge 
neutrality requires commensurate boson densities. Also, 
some disorder is unavoidably present, as proven by the 
analysis of the transition in finite magnetic fields in the 
same samples. Despite this disorder, in the vicinity of 
the zero field transition the temperature dependence of 
the resistivity and other quantities again do not indicate 
the presence of a Bose glass phase. 

In this letter we present a novel renormalization group 
analysis of the problem. We find that in dimensions 
d > 2 a direct insulator-superfluid transition takes place 
for weak disorder in these systems. Furthermore the dy- 
namical critical exponent z — d at the Bose Glass - Su- 
perfluid transition, where a statistical version of particle- 
hole symmetry develops. 

We consider the Hamiltonian: 



where Ui ( aj) annihilate (create) a boson at site i and 
fii — ajoi is the number operator, and fii = ^ + Ci, 
where the chemical potential /i controls the density of 
the bosons, and Ci is a random site energy, with distribu- 
tion P(e) over the finite support [—D, D]. 

For clean systems the phase diagram consist of a Su- 
perfluid phase (SF) and lobe-like Mott insulating phases 
(MI) with integer boson densities [g|. In the presence 
of disorder some site energies fall within the clean gap, 
shrinking the Mott lobes. It was proposed that a new, 
Bose Glass (BG) phase appears between the superfluid 
and the insulator, characterized by the absence of a SF 
order parameter, and a Mott gap. The excitations are 
presumed localized, with a finite density of states (DOS) 
at low energies and thus with a nonvanishing compress- 
ibility 1. 

We first present a mean field treatment of the prob- 
lem by chosing J,j = J/N, connecting all sites. In this 
formulation the asymptotics of the distribution of the 
site disorder plays a crucial role. It has recently been 
shown that for P(e) ~ (£) — |e|)", with a > 0, there is 
indeed a localization transition on the mean field level 
]l2[ . A physical motivation for such a choice is that spa- 
tial fluctuations of the field are likely to renormalize the 
customary uniform distribution to a smooth function in 
finite dimensions [Q. As this does not happen by it- 
self on the mean-field level, we instead build this feature 
in from the start by choosing a smooth unrenormalized 
distribution. We find that the critical behaviour at the 
MI-SF and the BG-SF transitions is independent of the 
details of the site-energy distribution. 

The evaluation of the free energy is achieved by per- 
forming a Hubbard- Stratonovich transformation on the 
kinetic term. In the N -^ oo limit the dynamic fluctua- 
tions are suppressed, thus the saddle point approximation 
becomes exact, yielding for the free energy 
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Siirn) = —^ifii + Ufii{ni — 1) — J{m*ai + ma\) , (3) 

where m takes its saddle-point value, m is proportional 
to the superfluid order parameter so it is zero in both 
the MI and BG phases. The MI phase is distinguished 
from the BG by the presence of a gap in its spectrum. 
The MI-BG transition is driven by the local collapse of 
this gap at some sites, which induces singular contribu- 
tions to /. While this is the dominant scenario at small 
J, for stronger kinetic couplings a direct MI-SF transi- 
tion is possible: long range fluctuations might generate a 
nonzero value for m before local instabilities would take 
place. In this case the free energy can be expanded in 
powers of m and a Landau-type action describes the MI- 
SF transition yielding the usual mean-field exponents. 

The advantage of employing a smooth distribution 
function of the disorder becomes apparent here. For the 
uniform distribution the MI-BG transition is absent for 
J > S, whereas it is manifestly present at some finite 
J for positive a's Q. 

We give credence to the above physical picture by de- 
veloping a new renormalization group (RG) scheme, sim- 
ilar to that of Ref. |1^ . As will be shown below, the soft 
core problem always scales to large interaction strengths, 
therefore we first present the procedure for hard core 
bosons and will discuss the effect of finite U later. Min- 
imizing the free energy yields the saddle-point equation 
at T^O 
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This sum is unfortunately plagued with singular terms 
with Hi ^ Q. The key concept of our RG technique is to 
avoid these singular terms by performing the sum step by 
step, integrating out those sites which have the highest 
excitation energies: E = maxj |/ii|. The contribution of 
these sites is then exactly absorbed into the renormaliza- 
tion of the parameters as: 

dm/dx = m(l — J/\/e^ + J'^m? 



(5) 



dJ/dx = j{z-l + .Jl\jE'^ + Pm? 

dD/dx = D(s- l/(a + l) 
dE/dx=^SE- D/{a + l) 



Here dx = — dlnN parametrizes the change in the num- 
ber of sites, and the dynamic critical exponent z re- 
lates this change to the rescaling the unit of energy SI: 
dJ7/Sl = zdN/N. These equations apply in the asymp- 
totic region, where P{e) ~ (Z? -I- e)". This region is 
reached by integrating out the high energy sites with 
Ei > in one step. Keeping the total bandwidth fixed, 
E — 1, requires z — D/{a + 1). 

The flow trajectories of the scaling equations are dis- 
played in Fig.l. for the case of m = 0. We find an at- 
tractive critical fixed point at D = J = 0. In the absence 



of disorder and hopping this point is naturally identified 
as the Mott Insulator. The critical fixed point a,t D = 0, 
J = 1 separates from the MI phase a region with run- 
away trajectories towards large J and vanishing disorder 
D. These trajectories are regularized by inducing finite 
values of m: thus they characterize a superfluid. This 
identification of a direct Mott-Superfluid transition for 
weak disorder is the central result of our paper. 
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FIG. 1. The renormalization group trajectories for a = 1. 

There is a critical point at Z) = 1, J = 0, which sends 
trajectories towards large D and J = 0, thus describing 
the MI-BG transition. As D increases, at the sites with 
the lowest local energy fii < 0, the cost of adding a par- 
ticle, \E — D\ eventually becomes equal to E. Since E 
is kept fixed at 1, this happens at Z? = 2, where half of 
the sites becomes empty. Subsequent RG steps should 
involve the elimination of the highest energy particle and 
hole excitations, from the two edges of the distribution 
simultaneously. From now on the trajectories remain in 
the D — 2 plane, since the width of the distribution 
(now determined by -E = D/2) is kept fixed. However 
the shape of the distribution of the site energies renor- 
malizes. Even if the original distribution is asymmetric, 
this procedure stretches the narrow section around the 
zero of energy up to the original bandwidth. Obviously 
the fixed point distribution is uniform, with half of the 
sites empty, representing a statistical particle-hole sym- 
metry, even though the original model did not possess 
this symmetry [|l4| . 

The fixed point distribution is thus characterized by 
a = 0, as well as E = 1, z = 1, yielding 
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The fixed point is at J = 0, I? = 2, m = 0, which 
is naturally identified as the BG-SF critical fixed point. 
Combining the equations for m and J we get dm/dJ ~ 
m(l/j2 - 1/J), yielding m ~ (1/J) exp(-l/J) [^. 
The following picture emerges for the flow trajecto- 



ries. For disorder values D > 1 the system first scales 
close to the MI-BG fixed point a,t D = 1, J = 0, which 
is characterized by the exponents, such as z = l/(ct + !)• 
But the flows continue towards the BG-SF fixed point 
at D = 2, J = 0. The corresponding critical behaviour 
becomes independent of a as the disorder distribution 
scales to a uniform function (a = 0). In particular, this 
restores the statistical particle-hole symmetry at the BG- 
SF fixed point, with z = 1. Here let us recall that our 
infinite-range hopping model is equivalent to the mean- 
field treatment of a finite-dimensional system. In a finite 
dimension d the dynamical exponent z relates the rescal- 
ing of the linear size, L, of the system to the energy-scale 
change f7'/f7 = {L' / ly . By writing [N'/N) = {L'/L^, 
from z = 1 we obtain z = d for the BG-SF transition. 
This result was first derived from the scaling of the com- 
pressibility in Ref. Q, but was recently debated [[l5[ . 
Finally, the special, multicritical fixed point at 13 = 1, 
J = a/{a + I) will be analyzed in a subsequent publica- 
tion. 
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FIG. 2. The phase boundary of the first Mott lobe for 



a = 1, for disorder values D 
represents the Bose-Glass. 
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Now we return to the case of soft core bosons. The 
novelty is that from the start particle and hole excita- 
tions at the same site have to be considered simultane- 
ously Formally this is implemented by expanding the 
saddle point equation in m, and then replacing l/|/ii| 
by [ni/Eh^i + [ui -\- 1)/Ep^i), where n.i is the number of 
bosons at site i in the ground state. Its value is deter- 
mined from 2U{ni — 1) < fii < 2Uni. The excitation 
energy of a hole is Eh,i — fn — 2U{ni — 1), and of a 
particle Ep,i — 2Uni — fit. Finally, U renormalizes, as 
dU/dx ~ zU . Let us observe that as long as the disor- 
der does not scale to zero, z = D/{a + 1) is positive, 
driving U to infinity. Thus for the finite D transitions 
the trajectories indeed scale to U = oo, thus we find the 
same fixed points, as in the hard core case. When the 
system scales towards small disorder, U increases, but 
saturates at a nonuniversal finite value. This stretches 



the hard core MI and the MI-SF fixed points into fixed 
lines. Therefore the fixed point structure of the soft core 
case is equivalent to that of the hard core case. 

We explicitly calculated the phase boundary of the 
Mott Insulator for soft core interactions by requiring that 
the coefficient of tti on the rhs. of the appropriately mod- 
ified saddle point equation equal 1. Since at the MI-BG 
phase boundary ^ is cither D or 2U — D, the straight 
separatrix at D = 1 in the D — J plane translates into 
straight sides of the Mott lobes, as shown in Fig. 2. It is 
worth noting that with increasing disorder the Mott lobe 
shrinks much faster along the /z axis, than along the J 
axis, in contrast to previous suggestions fs],^. 

Now we address the question of what happens on finite 
dimensional lattices. First the pure case will be reviewed, 
and then a small disorder turned on. We concentrate 
on the regions around the tip of the Mott lobes, where 
the transition happens at finite hopping J, so D can be 
regarded as a small parameter. 

In the pure case there are two types of transitions: exit- 
ing the MI by changing the chemical potential // (generic 
case) and by increasing the hopping J at the tip of the 
lobe (multicritical point) M. The generic transition is 
driven by particle or hole excitations. The dynamical 
exponent is z = 2 and the upper critical dimension is 
dc = 2. The gap linearly disappears with Sfi = fi — fic, 
requiring vz — 1, and i^ = 1/2 for any dimension d > 1. 
However the correlations inside the ground state do not 
diverge upon approaching the transition. On the other 
hand, the multicritical point possesses particle-hole sym- 
metry, setting the value of z to 1: the transition is of the 
type of a d -|- 1 dimensional XY model, hence dc = 3. 
The correlation length of the MI ground state does di- 
verge upon approaching the transition. 

Now we include site disorder D in the model. The 
relevancy of the disorder was extensively studied in past 
work. In classical models the perturbative Harris cri- 
terion is generally accepted to signal the irrelevancy of 
disorder |16|] . In their non-perturbative study, Chayes et 
al. proved that if the transition can be tuned by chang- 
ing the disorder around a finite disorder fixed point, then 
the correlation length exponent must obey v > 2/d [ [l7[ . 
This coincides with the Harris criterion for classical sys- 
tems, but is expected to apply for quantum transitions as 
well. In the quantum case, however, perturbative consid- 
erations may yield further bounds on exponents. Imag- 
ine for instance of tuning the system to criticality Here 
Lij{k) ~ fc^, by the definition of the dynamical critical 
exponent. In a finite size system this yields a level spac- 
ing for the first excited states of the clean Hamiltoninan 
5lj ^ L^^ . On the other hand, disorder introduces a mix- 
ing between these states, of the order of V^k' "^ L~'^/^. 
Second order perturbation theory shows that the disorder 
- induced shift of the energy levels is much less than the 
level spacing for d > 2z. Also, the perturbative change in 
the ground state wavefunction remains small, thus disor- 



der is irrelevant. All of the above arguments suggest that 
at the generic transition, where v = 1/2 and z = 2, weak 
disorder is irrelevant for d > dc — A. Thus the Bose-glass 
covers the sides of the Mott lobe only partially, leav- 
ing a finite region around the tip, where a direct MI-SF 
transition occurs, as in mean field. On the other hand, 
for d < 4 disorder becomes relevant, suggesting that the 
Bose-glass covers the sides of the Mott lobe completely. 

At the tip the Chayes criterion gives dc — 4 again, 
however the perturbative criterion yields dc — 2, since 
z = 1. These two results may be reconciled by consider- 
ing the following possibility. For D > one could argue 
that the randomness hinders ordering tendencies, thus 
the ordered superfluid phase should appear only at larger 
hoppings J. On the other hand, site - disorder is moving 
some energy levels inside the Mott gap. These sites may 
be able to support a superfluid even for hoppings smaller 
than Jc, the critical value for the clean system. Thus it 
is not inconceivable that these two competing influences 
cancel each other, and the phase boundary remains in- 
dependent of the disorder upto some finite D value. In 
this case the Chayes et al. bound does not apply, as its 
proof requires that the transition be driven by tuning the 
disorder. According to the other inequality then, the dis- 
order is relevant only for d < 2. Remarkably, all numeri- 
cal studies, performed at the tip of the Mott lobes, find 
a direct MI-SF transition MM- In the recent 2d work of 
Kisker and Ricgcr M it was found that for D < 0.4 J the v 
exponent at the tip of the lobe approximately equaled its 
clean value of I'adXY, violating the Chayes et al. bound. 
Correspondingly, the locus of the phase transition was 
consistent with the clean value, allowing for the possi- 
bility of systematic errors. Even more convincingly, the 
transition was characterized by genuinely new exponents 
around D « 0.4, such as z = 1.4, corresponding neither 
to the generic, nor to the XY criticality. Such a new 
multicritical behaviour is indeed expected, if a critical 
value of the disorder, Dc exists, such that for D > Dc 
the disorder is relevant. 

This scenario suggests that there can be a direct MI-SF 
transition even for finite disorder at the tip of the Mott 
lobe, for d > 2. One can rationalize this result by observ- 
ing that the correlation length diverges upon approaching 
the critical point from either side. One expects the dis- 
order to be screened and smoothed quite effectively by 
fluctuations of such large spatial extent. 

We wish to end on cautionary notes. Obviously these 
same numerical results can be in accord with the Chayes 
et al. bound, if the Bose glass phase is extremely slim, 
or if it manifests itself only on length scales, exceed- 
ing the biggest accesible system size. Such an anoma- 
lous behaviour however lacks a theoretical explanation 
to date. Additional subtleties associated with applying 
the Chayes et al. criterion around multicritical points 
were emphasized in Ref. 



of disorder on a lattice. We constructed the phase di- 
agram on the mean field level using an unconventional 
renormalization group scheme. At weak disorder a direct 
Mott Insulator- Superfluid transition takes place. We ar- 
gued that this transition is present for d > 4 in a finite 
region around the lobe-tips, whereas it survives down to 
two dimensions at the tip of the Mott lobes. Several nu- 
merical studies are consistent with this picture, as well as 
the limited experimental evidence on Josephson junction 
arrays. We proposed a possible reconciliation of these 
results with the well-known bound for v of Chayes et al. 
At strong disorder, we found that the Bose glass - Su- 
perfluid transition is characterized by statistical particle 
- hole symmetry and the exponent relation z = d on the 
mean field level. 
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